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This paper is devoted to the identification of the unknown smooth coefficient c 
PLh ' entering the hyperbolic equation c{x)dfu — Au = in a bounded smooth domain in M*^ 

. from partial (on part of the boundary) dynamic boundary measurements. In this paper 

I we prove that the knowledge of the partial Cauchy data for this class of hyperbolic PDE 

on any open subset T of the boundary determines explicitly the coefficient c provided 
that c is known outside a bounded domain. Then, through construction of appropriate 
test functions by a geometrical control method, we derive a formula for calculating 
the coefficient c from the knowledge of the difference between the local Dirichlet to 
^ ' Neumann maps. 
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■ 1 Introduction 

' In this paper, we present a differently method for multidimensional Coefficient Inverse 

Problems (CIPs) for a class of hyperbolic Partial Differential Equations (PDEs). In the 
k> i literature, the reader can find many key investigations in this kind of inverse problems, 

^ ; see, e.g. [2lll[6l[71[l6l[I71[22l|271l28llMll3l]and references cited there. L. Beilina and 

M.V. Klibanov have deeply studied this important problem in various recently works 
[His]. In [3], the authors have introduced a new globally convergent numerical method 
to solve a coefficient inverse problem associated to a hyperbolic PDE. The development 
of globally convergent numerical methods for multidimensional CIPs has started, as a 
first generation, from the developments found in \T8 \ \19 \ [20]. Else, A. G. Ramm and 
Rakesh have developed a general method for proving uniqueness theorems for multi- 
dimensional inverse problems. For the two dimensional case, Nachman [22] proved an 
uniqueness result for CIPs for some elliptic equation. Moreover, we find the works of 
L. Paivcirinta and V. Serov [231 ES] about the same issue, but for elliptic equations. 
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In other manner, the author Y. Chen has treated in [12] the Fourier transform of the 
hyperbohc equation similar to ours with the unknown coefficient c{x). Unhke this, we 
derive, using as weights particular background solutions constructed by a geometrical 
control method, asymptotic formulas in terms of the partial dynamic boundary mea- 
surements (Dirichlet-to-Neumann map) that are caused by the small perturbations. 
These asymptotic formulae yield the inverse Fourier transform of unknown coefficient. 

The ultimate objective of the work described in this paper is to determine, effec- 
tively, the unknown smooth coefficient c entering a class of hyperbolic equation in a 
bounded smooth domain in R*^ from partial (on part of the boundary) dynamic bound- 
ary measurements. The main difficulty appears in boundary measurements, is that the 
formulation of our boundary value problem involves unknown boundary values. This 
problem is well known in the study of the classical elliptic equations, where the char- 
acterization of the unknown Neumann boundary value in terms of the given Dirichlet 
datum is known as the Dirichlet-to-Neumann map. But, the problem of determin- 
ing the unknown boundary values also occurs in the study of hyperbolic equations 
formulated in a bounded domain. 

As our main result we develop, using as weights particular background solutions 
constructed by a geometrical control method, asymptotic formulas for appropriate av- 
eraging of the partial dynamic boundary measurements that are caused by the small 
perturbations of coefficient according to a parameter a. Assume that the coefficient is 
known outside a bounded domain 0, and suppose that we know explicitly the value 
of lim c{x) for x G fi. Then, the developed asymptotic formulae yield the inverse 



Fourier transform of the unknown part of this coefficient. 

In the subject of small volume perturbations from a known background material 
associated to the full time-dependent Maxwell's equations, we have derived asymptotic 
formulas to identify their locations and certain properties of their shapes from dynamic 
boundary measurements [13]. The present paper represents a different investigation of 
this line of work. 

As closely related stationary identification problems we refer the reader to [111 [T5| 
[22t [30] and references cited there. 

2 Problem formulation 

Let O be a bounded, smooth subdomain of M'^ with d < 3, (the assumption d < 3 is 
necessary in order to obtain the appropriate regularity for the solution using classical 
Sobolev embedding, see Brezis |9]). For simplicity we take Ofi to be C°°, but this con- 
dition could be considerably weakened. Let n = n{x) denote the outward unit normal 
vector to at a point on dQ. Let T > 0, xq € W^O, and let Q' be a smooth subdomain 
of Q. We denote by F CC as a measurable smooth open part of the boundary dQ. 
Throughout this paper we shall use quite standard L^— based Sobolev spaces to mea- 
sure regularity. 

As the forward problem, we consider the Cauchy problem for a hyperbolic PDF 




(1) 

(2) 
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where x(^) is the characteristic function of and ip E C°°{M. ) that ^(x) 7^ 0, Vx G fi. 

Equation ([T|) governs a wide range of apphcations, including e.g., propagation of 
acoustic and electromagnetic waves. 

We assume that the coefficient c{x) of equation ([T]) is such that 

\ca{x) = co(x) + aci(x) for x G r2, 
c(a;) = < , (3) 

[C2{x)= const. > for X G 

where Cj(x) G C^(r2) for i = 0, 1 with 

ci = in 17 \ iF, and M := sup{ci(x);x G n'}, (4) 

where Q' is a smooth subdomain of 0, and M is a positive constant. We also assume that 
a > 0, the order of magnitude of the small perturbations of coefficient, is sufficiently 
small that 

|ca(x)| > > 0, X G fi, (5) 

where c* is a positive constant. 

Suppose that the positive number C2 is given. In this paper we assume that the 
function c(x) is unknown in the domain ^l. Our purpose is the determination of c(x) 
for X G 0, assuming that the following function g{x,t) is known for the single source 
position xq G R'^\r2. Therefore, as done for the Dirichlet boundary conditions in [5], 
we set the Neumann boundary conditions: 

dv 

-Q:^\dQx{o,T) = 9{x,t). (6) 

The knowledge of c(x) outside of Q (c(x) = C2 in ]R'^\il), and the boundary function 
g{x,t) allow us to determine uniquely the function v{x,t) for x G as solution 

of the boundary value problem for equations ((l])-(2) with initial conditions in (2) and 
with the boundary conditions ([6]). Therefore, one can uniquely determine the function 
/(x,t) = v|anx(o,r)- 

Then, we can now consider an initial boundary value problem only in the domain 
X (0,T). Thus, the function v satisfying ([I])-(2) is en particular solution of the 
following initial boundary value problem 

{cad'^ - A)ua = in f7 X (0, T), 

Ua\t=o = V',dtUa\t=o = ip in f^, (7) 

Define u to be the solution of the hyperbolic equation in the homogeneous situation 
(a = 0). Thus, u satisfies 

{codf-A)u = mnx{0,T), 

u\t=o = (p,dtu\t=o = tp inO, (8) 
, u\dnx{o,T) = /• 
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Here G C°°{n) and / G C°°(0, T; C°°(5r^)) are subject to the compatibility conditions 
d^'f\t=o = (AV)bf7 and92'+V|t=o = (AV)lan, 1 = 1,2,... 

which give that dS]) has a unique solution in C°°{[0,T] x 17), see [H]- It is also well- 
known that ([7]) has a unique weak solution Ua G C'^(0, T; ^/'-'^(ri)) n C"^(0, T; see 

[21], [H]. Indeed, from [21] we have that -^len belongs to L2(0, T; ^^(aO)). 
Now, we define Fc := (917 \ T, and we introduce the trace space 

H^T) =|v G H^dn X (0,r)),t; = on x (0,r)}. 

It is known that the dual of H^{r) is H~^T). 

To introduce the local Dirichlet to Neumann map associated to our problem, we 
firstly define the function / = x(r)/ for {x,t) G dO, x (0, T), where x(r) is the 
characteristic function of T. Then, we have 

/ = /lrx(o,T) e^nd fem(T). (9) 

Therefore, we define the local Dirichlet to Neumann map associated to coefficient Ca 
by: 

A«(/) = ^|rfor/G#^(r), 

where Ua is the solution of ([7]). Let u denote the solution to the hyperbolic equation ([8]) 
with the Dirichlet boundary condition n = / on dO, x (0, T). Then, the local Dirichlet 

to Neumann map associated to cq is Ao(/) = T^lr for / G H^(T). 

on 

Our problem can be stated as follows: 

Inverse problem. Suppose that the smooth coefficient c(x) satisfies (3)- (4)- (5), 
where the positive number C2 is given. Assume that the function c(x) is unknown in 
the domain Q and / is given by Q. Is it possible to determine the coefficient Ca{x) 
from the knowledge of the difference between the local Dirichlet to Neumann maps 
Aq — Aq on r, if we know explicitly the value of lim (x) for x G ? 

To give a positive answer, we will develop an asymptotic expansions of an "ap- 
du 

propriate averaging" of — on T x (0, T), using particular background solutions as 

on 

weights. These particular solutions are constructed by a control method as it has been 
done in the original work [33] (see also [8], [TO], [21], [25] and [M])- It has been known 
for some time that the full knowledge of the (hyperbolic) Dirichlet to Neumann map 
{ua\dnx{o,T) '"^ ^^laQx(o,T)) Uniquely determines Conductivity, see [26], [31]. Our iden- 
tification procedure can be regarded as an important attempt to generalize the results 
of [26] and [31] in the case of partial knowledge (i.e., on only part of the boundary) of 
the Dirichlet to Neumann map to determine the coefficient of the hyperbolic equation 
considered above. The question of uniqueness of this inverse problem can be addressed 
positively via the method of Carleman estimates, see, e.g., [T7 1 [TO ] . 
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3 The Identification Procedure 



Before describing our identification procedure, let us introduce the following cutofi' 
function G C^{^) such that /3 = 1 on Q' and let r] G R*^. 

We will take in what follows (p{x) = e^^'^,ip{x) = —i\r]\e^'^'^ , and/(x,t) = e*'''^~*l''l* 
and assume that we are in possession of the boundary measurements of 

^ onrx(o,r). 

This particular choice of data 99, ip, and / implies that the background solution u of 
the wave equation ([5]) in the homogeneous background medium can be given explicitly. 

Suppose now that T and the part F of the boundary dil. are such that they geo- 
metrically control Q which roughly means that every geometrical optic ray, starting at 
any point x G at time t = hits T before time T at a non diffractive point, see [3]. 
It follows from [32] (see also [1]) that there exists (a unique) 5^ G Hq{0,T;TL'^{T)) 
(constructed by the Hilbert Uniqueness Method) such that the unique weak solution 

to the wave equation 

{cod^ - A)w^ = in Q X (0, T), 

dtWri\t=o = inil, (10) 

'"^»?lrx(0,T) = S'r?) 
, '^v\dQ\Tx(0,T) = ^' 

satisfies Wrj{T) = dtWr^iT) = 0. 

Let Orj G H^{0,T;L'^{T)) denote the unique solution of the Volterra equation of 
second kind 

rT 

dt9r,ix,t)+ / e''\^'^^^'-^\e.r,{x,s) -i\rj\dteri{x,s)) ds = g.r,{x,t) for X G r,t G (0,r), 

e^{x,0)=0 forxGT. 

(11) 

We can refer to the work of Yamamoto in [34J who conceived the idea of using such 
Volterra equation to apply the geometrical control for solving inverse source problems. 

The existence and uniqueness of this 0^ in H^{0,T; L^(T)) for any G M'^ can be 
established using the resolvent kernel. However, observing from differentiation of (llip 
with respect to t that 0^ is the unique solution of the ODE: 

9X - ^r, = e*l^l*at(e-^l''l*(7r,) forxGr,t G (0,r), . . 

9^{x,0) =0,dter,ix,T) = forxGT, ^ ' 

the function 0^ may be find (in practice) explicitly with variation of parameters and it 

also immediately follows from this observation that 6rj belongs to i?^(0, T; L^(r)). 

We introduce as the unique weak solution (obtained by transposition) in C^{0, T; L^(r2))n 
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C^{0,T;H ^(^2)) to the wave equation 



' {cod; - A)vrj = in 17 X (0, T), 
Vn\t=Q = in r2, 



dtVa,r,\t=o = iy ■ (??ci(x)e^''-") e L\n), 



(13) 

, '^»?l9nx(0,T) = 0. 
Then, the following holds. 

Proposition 3.1 Suppose that T and T geometrically control Q. For any rj G we 

have ^ 

I I gr,Ao{vrt) da{x)dt = / ci{x)e^''T'' dx. (14) 
Jo Jr Jw 

Here da{x) means an elementary surface for x G F. 

Proof. Let be the solution of ()13p . From [21] [Theorem 4.1, page 44] it follows that 

^oivr)) = "^~^lr ^ L^{^-,T]L^{T))- Then, multiply the equation {dj + A)t)^ = by 
and integrating by parts over (0, T) x 17, for any 77 G M*^ we have 



^ l^idf - A)v^wr, = ^ ^ V • i7]ciix)e'^-^)f3ix)e'^-^ dx - gr,^ 



Therefore 



since ci = on \ 17'. □ 

In term of the function Vr, as solution of (fTTj) , we introduce 

Ua{x,t) = u{x,t) +a'^ [ e~'^'^^''vr,{x,t-s)ds,xen,te{0,T). (16) 
Jo 

Moreover, for z{t) G C(f (]0,r[) and for any G Li(0, T; ^^(f^))^ ^e define 

v{x)= [ v{x,t)z{t) dt e L'^in). (17) 

JO 

The following lemma is useful to proof our main result. 

Lemma 3.1 Consider an arbitrary function c{x) satisfying condition ^ and assume 
that conditions and ^ hold. Let u, Ua he solutions of (0) and ^ respectively. 
Then, using \10^) the following estimates hold: 

\\Ua - '"||l°°(0,T;L2(Q)) < Coi, 

where C a positive constant. And, 

\\Ua — '«a||L°°(0,r;L2(n)) < C O.'^^^ , (18) 

where C is a positive constant. 
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Proof. Let ya be defined by 



^Va = Cadt{Ua — u) in Q. 



We have 



/ Cad^{ua-u)ya+ / V(Ua - u) • Vy„ = o / — Vu • V?/„. 
7n Jo Jo Co 

Since 



I V(n„ - u) ■ = - I Cadt{Ua - u){Ua - u) = -^-dt [ ( 

Jn Jn ^ Jn 



and 



/ Cadf{ua - u)ya = -^rdt [ |Vy, 



a I ) 



we obtain 



dt / |Vy„p + 9t / Ca{ua-uf = -2a / — Vu • Vy^ < C'a||Vya||L°°(o,r;L2(Q))- 
Jn Jn Jn co 

From the Gronwall Lemma it follows that 

\\Ua - w||l°=(0,T;L2{Q)) < Co. (19) 

As a consequence, by using (jl7p one can see that the function Ua — u solves the following 
boundary value problem 

{A{ua — u) = 0{a) in Q, 
{ua - u)\dn = 0. 
Integration by parts immediately gives, 

II grad {ua - ^i)||L2{n) = 0(a). (20) 

Taking into account that grad {ua — u)£ L°°(0, T; L^(r2)), we find by using the above 
estimate that 

||grad(u«-n)||i2(f2) =0(a)a.e. (0,r). (21) 
Under relation (I16p . one can define the function as solution of 

[ Aya = Cadt{ua - Ua) in Q. 
Integrating by parts immediately yields 

Jn Jn 

and 

/ V({ia - n„)Vya = -l/29t / Co,{Uo, - Uaf . 

Jn Jn 



To proceed with the proof of estimate (jlSp . we firstly remark that the function Uq 
given by (fTHl) is a solution of 

{cod^ - A)uo, = ia'^V ■ (r/ci(x)e^''-^)e-^l'?l* G L'^{n) inQx (0,r), 
Ua\t=o = v{x) in 17, 
dtUa\t=o = i^ix) in 17, 

„~, I ^iri-x—ilrilt 

Then, we deduce that Ua — Ua solves the following initial boundary value problem, 

{cad^-V- A)(u„-n„) = a'^V- (ci(x) grad ( / e"*l''l't;^(x, t - s) ds)) inl7x(0,r), 

Jo 

{Ua - Ua)\t=0 = inQ, 
dt{ua - Ua)\t=o = in 17, 

{Uq - Ua)\dnx{0,T) = 0- 

(22) 

Finally, we can use (f22]) to find by integrating by parts that 



JO JQ Jn 

which, from the Gronwall Lemma and by using ()2ip . yields 



I I ^ ' 

.Q, - Ua\\L^(o,T;L^{n)) S L/ O 



This achieves the proof. □ 



Now, we identify the function c{x) by using the difference between local Dirich- 
let to Neumann maps and the function 0^ as solution to the Volterra equation 
or equivalently the ODE ()12p . as a function of r]. Then, the following main result holds 



Theorem 3.1 Let rj G M'^, d = 2,3. Suppose that the smooth coefficient c{x) satisfies 
^M)-&-^M)■ Let Ua be the unique solution in C°(0, T; 7^^(17)) n C'^{0,T; L'^{n)) to the 
wave equation ([^ with (p{x) = e**^'^, ^l^{x) = —i\r]\e^'^'^ , and f{x,t) = e^'n-x-iHt ^ j^^i j 
he the function which satisfies Suppose that T and T geometrically control 17, then 
we have 

r I {er,+dte^df)iAa-Ao){f)ix,t) da{x)dt = a'^-^U^ / (c,-co)(x)e2*'?-" dx+0{a''+^) 

(23) 



a 



where Or/ is the unique solution to the ODE with defined as the boundary control 
in UOj) . The term 0{a'^~^^) is independent of the function ci. It depends only on the 
bound M . 
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~ du du 

Proof. Since the extension of (Aq, — Ao)(/)(x, t) to x (0, T) is — — ), then by 



, dury du , 



dn dn' 



conditions dt9r,{T) = and - ^)lt=o = 0, we have (A^ - Ao){f){x,t)\t=o = 0. 

Therefore the term 

r [ dtOrjdtiAa- Ao)if)ix,t) daix)dt 
Jo Jr 

may be simphfied as fohows 

' dtOr^dtiAa - Ao)(/)(x,t) da{x)dt 




Jr 

On the other hand, we have 




Jr 



d^9r,{Aa-Ao){f){x,t) da{x)dt. 

(24) 




Jr 



e^{K-K){f) + dt9r,dt{K-Ao){f) {x,t)da{x)dt = / / e^{AM) - Aa{Ua\r>ciO,T)))+ 



10 Jr L 

dt9rfdt{Aa{f) - Aa{ua\rx{o,T))) {x,t)da{x)dt+ 



, dv 



dvr 



e~'\'^\'^ix,t- s) ds + a'^dtOr^dt / e"*!''!" t - s) da{x)dt; 

on ' ^-^ 



dn 



where K{ua\rx{Q,T)) = Mf)+a'^ f e 'l''l'Ao(v^jr)(x, i - s) ds. 

JO 



Given that, O^j satisfies the Volterra equation (fT2]) and 



'0 



dtil e~'\'^\''^{x,t- s) ds) = at(-e-*l^l* S"^^'^{x,s) ds) 
./n an ./n an 



= i r? e 



i\ri\t I J\ri\s^]^ 



dVr 



[x,s) ds + -—^{x,t), 



dn ' dn 
we obtain by integrating by parts over (0, T) that 



'0 




'^{x,t){dter,+ r^,(s)e^l^l(*-^) ds) -i|r?|(e-^l''l*at0,(t)) f\^\^\^ ^ {x , s) ds) da{x)dt 



dn 



(x,t- s) ds + dte„dt I e-'^i^'^(x,t - s) ds]da(x)dt 
' dn J 

,dVr^ 
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Y dn 

-^{xM^A^ I (e^is) - i\r]\dterj{s))e'\^\^'-'Us) daix)dt 



dn 




10 Jr 

r-T 



^ ^ gr,{x,t)Ao{vr,\T){x,t) da{x)dt, 
and so, from Proposition 13.11 we obtain 

r [ k(A« - Ao)(/) + dtOr^dtiAa - Ao)(/)l ix,t)daix)dt 
Jo Jr^ J 



a%f [ ci{x)e^'T'' dx 
Jw 

i^(A„(/) - Aaiua\rx{o,T))) + dtO ^dt{Ao,{f ) - K{ua\rx(o,T))) da{x)dt + 0(a'^+^). 
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+ 



Jo Jr 



Thus, to prove Theorem 13.11 it suffices then to show that 



T 



9^{Aa{f)-K{ua\rx{o,T)))+dtOvdt{^M)-^a{ua\rx{o,T))) da{x)dt = 0(a'^+i). 



JT 

From definition (jl7p we have 

Ua-Ua= I {Ua - Ua)z{t) dt, 

Jo 

which gives by system (j22p that 

A{ua-iia) = [ Cadt{ua-Ua)z{t) dt+a'^ I V- (ci (x) grad ( / e"''^'^\''vn{x,t-s)ds))z{t) dt. 
Jo Jo Jo 

Thus, by (fT6|) and ([22]) again, we see that the function Ua — Ua is solution of 



-A{Ua - Ua) = - 
, {Ua - Ua)\dn = 0. 



Ca{ua — Ua)z"{t) dt + V ■ (ci(x) grad {ua — u)) in Q, 



(25) 

Taking into account estimate (llSp given by Lemma [3.H then by using standard elhptic 
regularity (see e.g. [H]) for the boundary value problem ([25]) we find that 



d 

\\-^{Ua - Ua)\\L'2{T) = 0{a 

~ ~ d 
The fact that Aa(/) - Aa{ua\rx{o,T)) ■= ^('"" ~ ^ L^^i^^T; L'^{T)), we deduce, 

as done in the proof of Lemma 13.11 that 



||A„(/) - A„(n„|rx(o,T))llL2(r) = 0{ 



a 



which implies that 



9^{Aa{f)-K{ua\rx{o,T)))+dtOr,dt{Aa{f)-Aa{ua\rx{o,T))) da{x)dt = 0(a°'+i). 
Jv^ J 

This completes the proof of our Theorem. □ 

We are now in position to describe our identification procedure which is based on 
Theorem 13.11 Let us neglect the asymptotically small remainder in the asymptotic 
formula ()23p . Then, it follows 

„—2ir]-x pT 



Ca{x)-Co{x) 



a' 



d-l 



\V\ 



Jr (^^+^*^''^*-)^^""^°^^-^~^^^'*^'^''^^^'^*'^'''^ ^ 



The method of reconstruction we propose here consists in sampling values of 
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at some discrete set of points r] and then calculating the corresponding inverse Fourier 
transform. 

In the following, it is not hard to prove the more convenient approximation in terms 
of the values of local Dirichlet-to-Neumann maps and Aq at /. 



Corollary 3.1 Let G M'^ and let f be defined by Suppose that T and T geomet- 
rically control 0,, then we have the following approximation 



Caix) « co(x)- 



a 



-2i-q-x pT 



Jr 



,^\v\tdt{e~'\^\'gr,{y,t)){A^-Ao){f){y,t)\da{y)dtdrj,x G Q, 

(26) 



where the boundary control is defined by U0\). 
f-T 



Proof. The term 

/o Jr 

to be interpreted as follows: 



dtOr]dt{Aa — Ao){f){x,t) da{x)dt, given in Theorem 13 -H has 



T 





dte^-dt{Ao,-Ao)if){x,t)da{x)dt 
dua du 



Jr 



d^e^ ■ {A^ - Ao)if){x,t)da{x)dt, 



because ^r?|i=T = and dt{- 



dn dn 



)\t=o = Q. In fact, in view of the ODE ([l2|), the 



term 



Jr 



Orj{Aa — Aq) + dfOrj ■ dt{Aa — Aq) f {x , t)da{x)dt may be simplified after 

integration by parts over (0, T) and use of the fact that 9^^ is the solution to the ODE 
(fT2]) to become 

- r [ e^\^\'dt{e~'\^\'gr,)-{A^-Ao){f){x,t)da{x)dt. 
Jo Jr 

Then, the desired approximation is established. □ 



4 Conclusion 

The use of approximate formula (j23p . including the difference between the local Dirich- 
let to Neumann maps, represents a promising approach to the dynamical identification 
and reconstruction of a coefficient which is unknown in a bounded domain(but it is 
known outside of this domain) for a class of hyperbolic PDE. We believe that this 
method will yield a suitable approximation to the dynamical identification of small 
conductivity ball (of the form z + aD) in a homogeneous medium in from the 
boundary measurements. We will present convenable numerical implementations for 
this investigation. This issue will be considered in a forthcoming work. 
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